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Abstract 

I present a general purpose Monte Carlo program for calculating the next-to-leading 
order corrections to arbitrary four-jet quantities in electron-positron annihilation. In 
the current version of the program, some subleading in color terms are neglected. As an 
example, I calculate the four-jet rate in the Durham scheme as well as the Bengtsson- 
Zerwas angular distribution at C(af ) and compare the results to existing data. 
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Four-jet events have been copiously produced at the Z-pole and at smaller energies 
at electron-positron colliders. These events allow the study of new aspects of QCD. For 
instance, the three-gluon coupling and the dependence on the number of light flavors, Nf, 
enter already at tree level. Therefore, this seems to be an ideal place for a measurement of the 
color factors |2|, |3[]. Even if there is no doubt about the correctness of QCD as the theory 
of the strong interactions, these measurements are not purely academic. In fact, the ongoing 
debate about the existence of light gluinos could be settled immediately by measuring Nf 
precisely enough (see e.g. ref. and references therein). Indeed, the addition of a massless 
gluino amounts to a change of Nf from 5 to 8. Also at LEP2 four-jet events will play an 
important role, since they are the main background for the W pair production. Thus, a 
next-to-leading order prediction for such quantities is highly desirable. 

In this talk I present some next-to-leading order results for the four jet fraction and 
some angular distributions. These results have been obtained with a program Q which 
can compute an arbitrary four-jet quantity at next-to-leading order. However, as will be 
discussed below, some subleading in color contributions have been neglected so far. 

Any next-to-leading order calculation involves basically two major steps. First, the corre- 
sponding one-loop amplitudes have to be computed and second, the phase-space integration 
has to be performed. This second step involves the cancellation of the real and virtual sin- 
gularities. I used the general version of the subtraction method as proposed in H to do 
the phase space integrals. Thus, no approximation at all has been made in this part of the 
calculation. 

The one-loop amplitudes which are needed for the calculation of four-jet events are 
e + e~ — ► qqq'q' and e + e~ — > qqgg. Recently, the amplitudes for four quarks in the final 
state have been computed |6|, [/J. However, the situation concerning the one- loop amplitude 
for the qqgg final state is less satisfactory. Unfortunately, only the leading in color terms 
are known so far ||. As a result, subleading in color pieces of the cross section can not yet 
be computed at next-to-leading order. Writing the color decomposition of any four-jet cross 
section as 



>) + El J b ) + El >) + El ff w + o (— 



I calculate a^ , ' b ' c ' < ^. 

Besides the one-loop amplitudes, the tree-level amplitudes for the processes e + e~ — > 
qqq'q 1 g and e + e~ qqggg are required for the computation of the real contributions. They 
have been obtained by several groups || [TOfl . In the program the results of ref are used. 

In the calculation of these amplitudes all quark and lepton masses have been set to 
zero. This is usually a very good approximation, although the 6-quark mass effects can 
yield considerable corrections O]. Unfortunately, the complete inclusion of mass effects at 
the order 0{oP s ) is presently out of reach, the main reason being the fact that the one-loop 
amplitudes are known only for the massless case. 

Besides the mass effects and subleading in color terms, three more contributions were 
neglected, although in principle their inclusion would not pose a problem. 

(1) Contributions coming from Pauli exchange. The corresponding 0{a 2 s ) terms are 
known to be numerically tiny and it is expected that the 0{a a s ) terms are numerically not 
very significant. 
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Figure 1: (a) Solid (dashed) lines represent the one- loop (tree-level) predictions for R4 in the 
Durham scheme for fi = y/s and a s = 0.118. (b) Solid (dashed) lines show the dependence of R4 
on the renormalization scale \x for the one-loop (tree-level) predictions in the Durham scheme, for 
y cut = 0.03. 



(2) Contributions proportional to the axial coupling a q of the Z° to quarks. Analogous 
terms have traditionally been omitted from 0{a 2 s ) programs, as they cancel precisely between 
up- and down-type quarks in the final state (for zero quark mass), and their contribution to 
the three-jet rate is at the percent level |fT2"| . 

(3) Contributions proportional to {Y. q v g) 2 y where v q is the quark vector coupling. These 
"light-by-glue scattering" terms do not appear at 0{a 2 s ) at all, have a partial cancellation 
from the sum over quark flavors, and contribute less than 1% to the 0{a z s ) term in the total 
cross-section [ p~3"f . 

I first present the results for the four-jet rate R4 = o"4_j e t/o"tot at next-to- leading order in 
a s in the Durham scheme [14]]. The y cn t dependence is shown in Fig. [I]. The solid (dashed) 
line represents the one-loop (tree-level) prediction. The renormalization scale \x has been 
chosen to be the center-of-mass energy y/s, the number of flavors Nf = 5 and a s = 0.118 [15 |. 
The data points are preliminary SLD data [R j and are corrected for hadronization. 

The truncation of the perturbative expansion for any physical quantity leads to a depen- 
dence of the theoretical prediction on the choice of the renormalization scale fi. The tree-level 
H dependence is much stronger for the four-jet rate than for the three-jet rate, because the 
former is proportional to a 2 instead of a s . As expected, this strong renormalization-scale 
dependence is reduced by the inclusion of the next-to-leading order contribution. Fig. p] also 
plots the //-dependence of _R 4 at tree-level and at one-loop for y cut = 0.03. 

Four-jet angular distributions [[HJ have been measured by several collaborations Q, |3[ . 
The general procedure is to choose a certain jet definition. Then, in the case of a four- 
jet event, the jets are ordered according to their energies such that E 1 > E 2 > E 3 > £4. 
Usually, the most energetic jets are associated with the primary quarks whereas the remaining 
two jets either origin form a quark or gluon pair (at tree level). This can be exploited to 
construct angular variables which have a completely different shape for the four-quark and 
the two-quark-two-gluon final state. Since the two final states are proportional to different 
color structures one can attempt to measure the various color factors and in particular the 
number of light flavors Nf. Unfortunately, the four-quark final state is strongly suppressed. 
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Figure 2: Bengtsson-Zerwas distribution at tree level (dotted), one-loop (solid) and one-loop 
with Nf = 8 (dashed) compared to (a) OPAL and (b) ALEPH || data (histograms), which are 
corrected for detector and hadronization effects. 



As a result, the full distributions are not very sensitive to Nf and the error on the measured 
value of Nf is accordingly large. 

An advantage of the angular distributions lies in the fact that one does not need to worry 
about large logarithms coming from a particular choice of the renormalization scale. The 
reason for this can easily be understood. At tree level, the strong coupling constant appears 
only in an overall prefactor. Since these distributions are normalized, the value of a s and 
thus the choice of the renormalization scale \i has no influence at all on the result. Only the 
inclusion of the one-loop corrections introduces a extremely mild /i-dependence. 

As an example, I consider the Bengtsson-Zerwas angle, Xbz, and compare the next-to- 
leading order prediction with the two most recent analyses of OPAL and ALEPH 0. |3|. 
In jets were defined according to the JADE scheme with y cut = 0.03, whereas in || the 
Durham jet algorithm with the E0 recombination scheme was used and y cut was chosen to be 
0.008. Note also that the two experiments use different normalizations. Fig. compares the 
tree-level (dotted) and next-to-leading order (solid) predictions to the data which have been 
corrected for detector and hadronization effects. The dotted line can hardly be seen because 
it nearly coincides with the solid line. The theoretical curves have been obtained by binning 
Xbz into twenty bins. This rather fine binning results in a somewhat larger statistical error, 
which is of the order of 2% for the shown curves. In order to illustrate the mild dependence 
on Nf I plotted also the one- loop results for Nf = 8 (dashed). Although this dependence 
and thus the precision on the measurement of Nf may be enhanced by additional cuts or by 
6-quark tagging |18" |, Fig. |2] shows that it is very difficult to get a precise measurement of Nf 



from angular distributions alone. 
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